Universal mesoscopic statistics and the localization of light 
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We follow the evolution with sample thickness, of intensity statistics for localized light transmitted 
through layered media in a crossover from one to three dimensions occasioned by transverse disorder. 
The probability distribution of intensity changes from one dimensional to a mixture of a mesoscopic 
function of a single parameter, the "statistical conductance," and a distribution of intensity for 
Gaussian waves. This suggests that the change to a universal statistics beyond ID is associated 
with the topological change in the spatial field distribution. 
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The flow of energy in random systems depends upon 
the closeness to the threshold for Anderson localization 
separating diffusing and localized waves. [l|, [if On one 
side of this divide, average transport may be described 
by the diffusion of particles, such as electrons or photons. 
On the other side, propagation is strongly suppressed 
and proceeds via the excitation of spatially localized reso- 
nances of the underlying quantum mechanical or classical 
waves, d, 0, m The nature of transport in an ensemble of 
nondissipative random samples can be gleaned from the 
scaling of the ensemble average of the dimensionless con- 
ductance, g = G/(e^//i), which depends only upon the 
average value of g and the dimensionality of the scaling. 
^] Following convention, g and its ensemble average will 
be denoted by the same symbol. Localization is achieved 
for ^ < 1. The Landauer relation shows that g is equiv- 
alent to the optical transmittance, g = T = ^^^^Tab^ 
where Tab is the transmission coefficient for an incident 
transverse mode a into the outgoing mode or a point 
on the output surface, b. [6] Thus localization may be 
achieved in principle for light as well for electrons. Lo- 
calizing classical waves is of particular interest because 
it is possible to study both local and global fluctuations 
of transmission in random ensembles and to study An- 
derson localization without the complication of electron 
interactions. 

Photon localization has been demonstrated in one [7] 
and two dimensions, [8, 9] as well as in quasi- ID samples, 
(iq| in which the length of a locally three dimensional 
sample with reflecting side walls greatly exceeds the di- 
ameter. However, it is difficult to achieve and demon- 
strate definitively in random three dimensional samples 
without structural correlation. ^ In contrast to elec- 
trons, which exhibit 5-wave scattering, the leading order 
of scattering for light is p-wave so that the scattering 
cross section is at a maximum when the scale of the scat- 
tering elements is of the order of the wavelength of light, 
D ^ X. This limits the density of scatterers and con- 
sequently the minimum achievable transport mean free 
path for photons, £. As a result, scattering may not be 
strong enough to satisfy the loffe-Regel criterion for lo- 



calization in three dimensions, ki = 1^ where k = 2tt/\ is 
the wave vector. [12] At the same time, the exponential 
scaling of optical transmission [13] is not an unambiguous 
sign of localization since absorption in diffusive samples 
also leads to the exponential falloff of transmission, [ij] 
Because of large fiuctuations in mesoscopic systems, 
which increase dramatically in the localization transi- 
tion, it is essential to measure the statistics of transport. 
The probability distributions of total transmission and 
intensity normalized to their respective ensemble aver- 
ages, P{sa = Ta/{Ta)) and P{s 

ab — Tdb /{Tab)), are di- 
rectly related in quasi- ID systems because the transmit- 
ted wave is thoroughly mixed so that intensity statistics 
are uniform across the sample output. These distribution 
were calculated in the diffusive limit for nondissipative 
saniples and expressed in terms of a single parameter, 
g. [H, [l^ Since g = 2/3var(5o) in this limit, the distri- 
bution of P{sa) depended only upon the variance of the 
distribution. The functional forms of P{sa) and P{sab) 
were found to be valid even in strongly scattering and 
dissipative quasi- ID samples once the "statistical con- 
ductance", g^ = 2/3var(5a) is substituted for g. [l7| g^ 
scaled inversely with L for g^ > 1 and fell exponentially 
for < 1, [10] as was originally predicted for the local- 
ization parameter g. [2] 

Recently, the intensity distribution was compared to 
measurements of pulsed and steady ultrasound in a slab 
of randomly positioned aluminum beads sintered to form 
an elastic network. [18] Measurements of P{sab) were 
well fit by the distribution in Eq. ([T]) below at different 
frequencies with a minimum in g' 0.8. Agreement with 
the distribution found in quasi- ID samples is intriguing 
in this 3D system since the probability distribution of in- 
tegrated transmission depends crucially on the area over 
which intensity is averaged, but there is no specific area 
over which Sa or g' can be unambiguously defined in the 
slab geometry. 

Since dimensionality is crucial to scaling, [2] a study 
of transmission statistics in a continuous change of di- 
mension may elucidate the general character of statistics 
in random media. A change of dimensionality can be 
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induced in random layered media with transverse dis- 
order. [19] Measurements of average transmission in 
stacks of glass cover slips showed a change from one- 
dimensional scaling, with an asymptotic exponential scal- 
ing for (T(L)) ~ exp(— where ^ is the average 
localization length of intensity, to an inverse scaling of 
transmission, characteristic of diffusion. Disordered lay- 
ered systems are ubiquitous and present opportunities 
and challenges in diverse settings such as multilayered 
Bragg gratings, acoustic scattering in the earths crust, 
and anomalous conductivity anisotropics [2q| . 

In this Letter, we follow the evolution of P{sab) for 
plane wave excitation for localized light in a crossover 
from one to three dimensions in layered media induced by 
transverse disorder. In thin samples in which light does 
not spread beyond a single coherence area, P{sab) is in 
accord with ID simulations, and falls sharply to zero as 
Sab and for Tab > 1. In thicker samples, the spatial 
intensity pattern undergoes a topological transformation 
in which phase singularities are formed at intensity nulls 
as the lateral spread of the wave exceeds the field corre- 
lation length. P{sab) changes to the distribution given 
in Eq. (p!a|) below, which was originally found for quasi- 
ID system. This distribution can be written formally 
as a mixture of a mesoscopic distribution of Sa, P{sa)^ 
and a distribution of intensity for Gaussian waves, even 
though there is no area for which the distribution of in- 
tegrated intensity is given by the quasi- ID expression for 
P{sa)' These measurements show that mesoscopic in- 
tensity statistics have a universal form beyond ID and 
provide a definitive test for Anderson localization. 

The sample is composed of glass cover slips with inter- 
vening air gaps. The cover slips are 22 x 22 mm^ with av- 
erage thickness of 150 /im and standard deviation of the 
wedge angle of ~ 0.05°. The refractive index of glass is 
1.522 at the 632.8 nm wavelength of the polarized single- 
frequency helium- neon laser used in the experiment. The 
first surface of the sample is at the waist of the normally- 
incident 5-mm-diameter beam. A tapered optical fiber 
with mode field diameter of 4 jam was centered in the 
beam and placed just behind the stack. The signal was 
detected by a PMT. The sample was translated in the 
layer plane over a 2 x 2 — mm^ area on a 10 — fim grid to 
yield a speckle pattern equivalent to that for an incident 
plane wave. For each sample thickness, speckle patterns 
were obtained for 10 different arrangements of cover slips. 
Examples of speckle for 20 and 60 cover slips are shown 
in Fig. 1. 

Measurements of var(sa5) versus number of cover slips 
L are shown in Fig. 2 and compared to ID transfer ma- 
trix simulations. The measured variance is in excellent 
accord with ID simulations in thin samples, but depar- 
tures can be seen for L > 30. Beyond this thickness, the 
average transverse spread of the wave, a^, exceeds the 
field correlation length, 1/Sk±^ where 5k± is the width of 
the distribution of k±. [19] Thus, for 5k±(j± > 1, wave 




FIG. 1: (Color Online) Typical speckle patterns measured at 
the output of the samples with (a) L = 20 and (b) L = 60. 
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FIG. 2: (Color Online) Comparison of measured var(sab) 
changing with L to ID-simulation. 



propagation is no longer one-dimensional and transmis- 
sion is not determined solely by the longitudinal variation 
of the index of refraction along a line through the sam- 
ple but also refiects lateral scattering due to transverse 
inhomogeneity. The measured values of var(5a6) reach a 
maximum value of 3.25 at L ~ 80 while the simulated 
results increase exponentially. This may be compared to 
the value at the localization threshold in quasi- ID deter- 
mined by the condition, a' = 1, var(sa6) = 2var(5a) + l = 
4/3^' + 1 = 7/3. [l5|, 17|. The question arises as to 
whether the value = 0.63 inferred from the above re- 
lation with YdiT{sab) = 3.25 signifies that the wave is lo- 
calized in the slab geometry of a layered medium. This 
question may be addressed by examining the full inten- 
sity distribution as L is changed and comparing the re- 
sults to the quasi- ID distribution. 

The relationship between the statistics of Sab and Sa 
in quasi- ID follows from the statistical independence of 
= Sab/sa = and Sa and the Gaussian field statistics 
in single configurations. pj|, [IJ] Fluctuations of inten- 
sity in quasi- ID samples may therefore be represented as 
a mixture of P{sa) and P{s'^-^). Since P{s^ab) ^ nega- 
tive exponential distribution for polarized transmission. 



P{s'^^) = exp(— 5^^), the intensity distribution is given 
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by, 



P{Sab) = [ —exp(-—]P{Sa)dSa, 
Jo \ Sa / 

^ nioo 
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=exp -g'W[^\^q^^] . 
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The distribution P[^Sa) is found from RMT [l5| and di- 
agrammatic [16] calculations in the weak scattering limit 
in the absence of dissipation to be a function of a single 
parameter g. The distribution has been found to ap- 
ply more broadly in quasi- ID to dissipative and strongly 
scattering systems once g' = 2/3var(5a6) is substituted 
for g in the expression for P(sa), [17] 

Measurements of P{sab) in layered samples are shown 
for L = 20, 30, 60 and 100 in Fig. 3 and compared to 
ID simulations and to Eq. ([Taj). P{sab) for L = 20 is in 
excellent agreement with ID simulations up to Sab = 2.5. 
The maximum value of Sab^ above which ID simulations 
of P{sab) vanishes, which is 3.4 for L = 20, corresponds to 
full transmission. Tab = 1- Higher values of transmission 
occur, however, because of the transverse spread of the 
beam due to nonuniformity in the layers, which does not 
exist in ID. P{sab) is expected to be high for small values 
of Sab for localized waves in ID samples since transmis- 
sion is low when the light is off resonances with localized 
modes and when modes are peaked away from the cen- 
ter of the sample. However, as seen in the inset in Fig. 
3(a), P{sab) drops precipitously towards zero as Sab 0. 
The absence of nulls in transmission is a distinctive fea- 
ture of ID media and follows from the finite transmission 
coefficient for flux through the i/i -\- 1 interface between 
parallel layers, 2nini-^i/{ni + 71^+1)^ . This requires that 
the component of flux propagating towards the sample 
output does not vanish in a layer if it is nonzero in the 
preceding layer. Nulls in intensity of transmitted light 
can occur, however, in layered samples once the distri- 
bution of transverse wave vectors, k±^ is broadened by 
scattering from layers with nonpar allel surfaces. 

For L = 30, the intensity distribution falls between the 
ID and quasi- ID distribution, whereas P{sab) for L = 60 
and 100, are in agreement with Eq. ([Taj) using the value 
of g' obtained from the relation, g' = [var(5a6) — l]/2, 
within the limits set by the mode field diameter of the ta- 
pered optical pickup fiber. The agreement with Eq. ([Taj) 
in this highly anisotropic medium in a three-dimensional 
slab geometry is surprising since the statistics of propaga- 
tion in quasi- ID flow from the complete mixing of waves 
at the sample output which gives uniform transmission 
statistics across the output. But, in the slab geometry, 
the mixture of transmission channels as well as the quasi- 
modes through which the wave is coupled to the output 
varies along the output. Further, whereas the statistics of 
Sa depend critically on the sample cross section in quasi- 
ID samples, there is no well-defined area over which Sa 
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FIG. 3: (Color Online) Comparison of the measured P{sab) 
to ID-simulation and Eq. ([Ta)) with same var(sab)- Samples 
studied are with L = 20 for (a), 30 for (b), 60 for (c) and 100 
for (d). Inset of (a) shows the details of P{sab) for small Sab, 
demonstrating the absence of phase singularities at L = 20. 



can be defined in the slab geometry. The distributions 
of integrated transmission through a number of different 
circular areas at the output of the sample with L = 60 
are shown in Fig. 4 and compared to the distribution of 
P{sa) for the value of g^ inferred from the measurements 
of P{sab)' None of the curves for integrated transmission 
correspond to the distribution predicted by Eq. ([Taj) for 
P{sa)- Yet, Eq. ([Taj) quantitatively describes intensity 
statistics, once Sa is treated as a variable of integration 
rather than as a random variable. This suggests that g' 
is a generalized measure of mesoscopic fluctuations not 
directly related to the statistics of integrated transmis- 
sion. It may be related to the probability of crossing of 
Feynman paths at scattering centers within the medium 
or the number of electromagnetic quasimodes within the 
medium contributing to the field at the output. We find 
that for L > 40, ^Ms equal to the product of the number 
of coherence areas in an area over which the incident field 
spreads at the output, N = {Sk±(j±)'^ , and the ensem- 
ble average transmission coefficient, g' = N{Ta) , where, 
(Ta) = N{Tab). Since, g = Eab^ab = N{Ta) , in this 
case, g' = g. 

The change of P{sab) with increasing thickness is a 
consequence of the change of the spatial distribution of 
intensity. Beyond ID, P{sab) is obtained by mixing the 
mesoscopic distribution, P(sa), with a distribution of in- 
tensity for a random Gaussian process [Eq. ([Taj) ]. The 
change of the spatial distribution of transmitted flux from 
ID to higher dimension is a distinctive topological trans- 
formation accompanied by the formation of a speckle pat- 
tern with intensity nulls which lie at the center of vortices 
of electromagnetic flux. [22I, 113] Such optical vortices 
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FIG. 4: (Color Online) Probability distributions of intensity 
integrated over circular areas with different radii for L = 60. 
The results are compared with P(s a ) given by Eq. (Ilbl) . 

form a network which underpins random Gaussian fields, 
[m Low intensity values within the speckle pattern pre- 
dicted by Eq. (p!a|) are associated with points within the 
vortex core. Gaussian statistics arises from the sum of 
numerous statistically independent contributions to the 
output field. As the number of transverse modes con- 
tributing to the field at a point increases, the assumption 
of underlying Gaussian statistics for intensity are more 
nearly satisfied. 

In conclusion, we have examined the statistics of light 
localized longitudinally by multiple scattering and lat- 
erally by variations in the transverse and longitudinal 
disorder in the sample which create a pattern of trans- 
mission resonances. Localization in this anisotropic 3D 
sample is due to the restricted phase space for k which 
is a narrow cone around the forward and backward di- 
rections. The statistics of light transmitted through lay- 
ered media with transverse disorder is ID as long as the 
transverse spread of the light is smaller than the coher- 
ence length. The character of intensity statistics changes, 
however, once phase singularities appear in the transmit- 
ted field. Intensity statistics beyond ID are then a mix- 
ture of the intensity distribution for a Gaussian field and 
a mesoscopic distribution, which depends only upon the 
mesoscopic variable g\ The similarity of local and global 
statistics of intensity and evolution of the speckle pat- 
tern in quasi- ID [21] suggests that a universal statistics 
should also exist for aspects of speckle change for any 
dimension beyond ID. Thus, the statistics of generalized 
velocity of phase singularities with frequency shift or in- 
ternal change of the sample, for example, should be a 
mixture of the velocity distribution for a Gaussian field 
and a mesoscopic distribution in terms of a parameter 



proportional to g' . The excellent agreement of Eq. ([T]) 
with measured intensity distributions even in the case 
that the speckle pattern is not generic and reflects the 
underlying disorder, as may be seen in Fig. 1, suggests 
that Eq. ([T]) is robust and holds for any form of disorder 
above ID. 
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